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ABSTRACT 



. Context. Radiative transfer calculations have predicted intensity enhancements for optically thick emission lines, as 

opposed to the normal intensity reductions, for astrophysical plasmas under certain conditions. In particular, the results 
, are predicted to be dependent both on the geometry of the emitting plasma and the orientation of the observer. Hence 

I in principle the detection of intensity enhancement may provide a way of determining the geometry of an unresolved 

, astronomical source. 

Aims. To investigate such enhancements we have analysed a sample of active late-type stars observed in the far ultra- 
violet spectral region. 

Methods. Emission lines of O vi in the FUSE satellite spectra of eEri, II Peg and Prox Cen were searched for intensity 
^ , enhancements due to opacity. 

■ Results. We have found strong evidence for line intensity enhancements due to opacity during active or flare-like ac- 

' tivity for all three stars. The O vi 1032/1038 line intensity ratios, predicted to have a value of 2.0 in the optically thin 

I case, are found to be up to ~ 30% larger during several orbital phases. 

Conclusions. Our measurements, combined with radiative transfer models, allow us to constrain both the geometry 
of the O VI emitting regions in our stellar sources and the orientation of the observer. A spherical emitting plasma can 
be ruled out, as this would lead to no intensity enhancement. In addition, the theory tells us that the line-of-sight to 
^ ' the plasma must be close to perpendicular to its surface, as observations at small angles to the surface lead to either 

m : no intensity enhancement or the usual line intensity decrease over the optically thin value. For the future, we outline 

a laboratory experiment, that could be undertaken with current facilities, which would provide an unequivocal test of 
'/^ ' predictions of line intensity enhancement due to opacity, in particular the dependence on plasma geometry. 



Key words, stars: activity — stars: late-type — opacity 



1 . Introduction impact excitation rates, which is 2.0 (Aggarwal & Keenan 

2004). Hence deviations of the observed line ratio from 

Deriving the physical extent and geometry of spatially un- "P^i'^^^y ^^^^^ "Pti'^^l ^^^P* ^"^^ h™'^^ 

resolved stellar corona and transition regions (TRs) is an Parameters of the emitting plasma to be derived. Using 



H ; important undertaking. Previous studies seeking to derive "i^*"'^' ^^^^''^1 ^^^^11^^ ^ave found TR scale heights 

. scale heights for stellar TRs (Mathioudakis et al. 1999; ^^1^^^" \0-100 km for active stars including AUMic 

Bloomfield et al. 2002; Christian et al. 2004) used the ^^1°°™,^!'^ 2002), the flare star Prox Cen (Christian 

relation of optical depth to physical parameters, includ- 2004) and AD Leo (Christian et al. 2006). 

ing electron density and pathlength. Transitions such as We note that 1032/1038 intensity ratios of Ovi with 

1032 A (2s ^S-2p ^P3/2) and 1038 A (2s ^S-2p ^Pi/2) values ranging from 1-4 have been measured in so- 

of O VI have a common lower level, and under optically lar spectra obtained by the Ultraviolet Coronagraph 

thin conditions in a collisionally-excited plasma their ex- Spectrometer (UVCS) on the SOHO satellite (see, for 

pected line intensity ratio scales as that of their electron example, Nakagawa 2008). However, these observations 
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were obtained at large heights above the solar sur- 
face, where the electron densities are very low (around 
10^ cm"^ or less; Ko et al. 2006), and the O vi line emis- 
sion has a significant component from the resonant scat- 
tering of chromospheric O vi radiation and/or the absorp- 
tion and subsequent re-emission of Doppler- shifted C n 
1036.3 and 1037.0 A photons by Ovi 1038 A (Kohl & 
Withbroe 1982). The O vi line intensities decrease rapidly 
with increasing height above the solar surface, with for 
example a factor of 10 reduction between 1 and 1.1 Rq 
(Nagakawa 2008). Hence for an unresolved stellar source, 
such as those considered here, the O vi line emission will 
be dominated by the high density coUisionally-excited 
component and will have an optically thin 1032/1038 in- 
tensity ratio of 2.0. 

Recently, theoretical work has indicated that in cer- 
tain instances the intensity ratio of an optically thick to 
optically thin line in a collisonally-excited plasma could 
increase (rather than the expected decrease), due to the 
fact that an ion in the upper state of the transition can 
be pumped in the optically thick case by photons travers- 
ing the plasma at many different angles (Kerr et al. 2004, 
2005). In addition, these studies predicted that the degree 
of line enhancement for the optically thick transition de- 
pends on both the geometry of the emitting plasma and 
the orientation of the observer. This is an important result, 
as it potentially provides a way of determining the geom- 
etry of a spatially unresolved plasma using purely spec- 
troscopic means. Some evidence for such a line enhance- 
ment was recently found for the active star EVLac us- 
ing observations of the 15.01 and 16.78 A lines of Fe xvn 
from the XMM-Newton satellite. Rose et al. (2008) mea- 
sured a 15.01/16.78 intensity ratio of 2.50 + 0.25 (Icr er- 
ror), compared to a theoretical optically thin value of be- 
tween 1 .75-1 .93, indicating some enhancement in the op- 
tically thick 15.01 A transition. However, we should note 
that there are uncertainties in the atomic data for Fe xvn 
and the theoretical fine ratios are also sensitive to the 
adopted plasma conditions (see GiUaspy et al. 2011 and 
references therein). 

In the present work, we extend our study to 
far-ultraviolet observations of active late-type stars. 
Specifically, we use moderate to high resolution spectra 
from the Far Ultraviolet Spectroscopic Explorer (FUSE) 
satelhte to search for line intensity enhancements due 
to opacity in the atmospheres of several late-type stars. 
We examine the 1032 and 1038 A emission lines of Ovi, 
where the transitions have a common lower level and the 
theoretical value of the optically thin 1032/1038 ratio in- 
volves simpler and more reliable atomic physics, plus has 
a lesser dependence on the plasma parameters, than is the 
case for Fe xvii. In Section 2 we describe the FUSE ob- 
servations and data reduction techniques, while Section 3 
provides details of the theory used in the analysis. Our re- 
sults are presented in Section 4, including measurements 



of the emission Une fluxes, and the determination of den- 
sities and pathlengths for the emitting plasmas under con- 
sideration. In Section 5 we discuss our line ratios and im- 
plications for the geometry of the emitting region, and 
also outline some possibiUties for future work, including 
a laboratory experiment to test the theory behind this re- 
search. 



2. Observations 

The FUSE satellite operated successfully from its launch 
in June 1999 until October 2007. It consisted of four 
co-aUgned prime-focus telescopes, with two having SiC 
coatings and optimised for the 905-1 105 A region, and 
the others using LiF coatings to cover 987-1 187 A. 
Holographically-ruled gratings were used to disperse 
fight onto two separate microchannel plate detectors 
(MCPs), each with two independent segments. The spec- 
tral resolution for the HIRS aperture was ~ 20,000 and 
~ 12,000-15,000 for LWRS. Details of the FUSE instin- 
ment and in-orbit performance may be found in Moos et 
al. (2000) and Sahnow et al. (2000a,b). 

2.1. Sample selection 

The FUSE observations and datasets were obtained from 
the Multimission Archive at Space Telescope (MAST). 
Quick-look spectra for the Ovi lines (1032 and 1038 A) 
were investigated and stars with possible line ratio en- 
hancements selected. These include the active star eEri 
and the RS CVn II Peg. In addition, we employ Ovi line 
ratio results for Prox Cen from Christian et al. (2004). 
Table 1 summarises the observational datasets for the stel- 
lar sample used in the current study. 



2.2. Data analysis 

The raw FUSE data files were re-processed with the lat- 
est version of the FUSE caUbration pipeUne (currently 
CalFUSE V3.2). Spectra were extracted from the appro- 
priate aperture (HIRS or LWRS) and background sub- 
tracted, flat-fielded, wavelength and flux calibrated in the 
standard CalFUSE reduction (Dixon et al. 2007). The 
eEri spectra were obtained in the high resolution mode 
(HIRS), and II Peg in the LWRS mode. 

Spectra were extracted as a function of intensity and 
phase using the new IDF.CUT routine. The first 25 ksec 
of the eEri observations were a factor of 25 brighter than 
the remaining ones (see Fig. 1), which had a median count 
rate of 0.009 counts/sec. This initial 25 ksec was divided 
into 8 phase bins. Most of the II Peg observations were 
active and this entire dataset was divided into 1 1 phase 
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Source 


Program ID 


Exp" 


Date" 


aperture 


Spectral type 


Comment 


Prox Cen 


D122 


45.4 


2003-04-05 


LWRS 


M5.5 Ve 


Flare star 


eEri 


C165 


73 


2003-12-21 


HIRS 


K2 V 




II Peg 


P179 


58 


2000-11-23 


LWRS 


KOIV 


RSCVn 



Total exposure time in ksec. 

'Observation date (earliest date given if more than 1 sequence). 



bins. Extracted spectra were fitted using IRAF and DIPSO 
routines. 

As in our previous work (Bloomfield et al. 2002; 
Christian et al. 2004, 2006) we measured emission line 
fluxes by fitting Gaussians to the observed line profiles. 
In general, a single Gaussian profile was sufficient to fit 
both the O vi transitions. 

Light curves extracted in the lines of Ovi (1032 + 
1038 A) are shown in Figs. 1 and 2 for eEri and II Peg, 
respectively. 

3. Theory 

In this section we extend the previous theoretical work 
by Kerr et al. (2005) on the variation of line intensity 
enhancement with optical depth, and its relation to the 
plasma geometry. The research of Kerr et al. made sev- 
eral assumptions, including uniform plasma conditions 
and specific geometries. However, the lines they consid- 
ered involved different kinetic pathways to the two transi- 
tions discussed here and so we need to alter their analysis 
in this respect. Specifically, Kerr et al. investigated one 
optically thick and one optically thin line which did not 
involve any common states. In our case, the two lines of 
Li-like O vi at 1032 and 1038 A involve a common lower 
level (the ground state 2s ^S), with oscillator strengths and 
hence optical depths differing only by a factor of two. The 
densities and temperatures involved are such that the ki- 
netic system can be considered to be in coronal steady- 
state, with the upper two levels (2p ^Pi/2 and ^Pa/i) not 
connected kinetically. If an analysis similar to that under- 
taken by Kerr et al. is applied to this problem, and not- 
ing that the electron collisional excitation rates for the 
two transitions 2s ^S-2p ^¥3/2 (1032A) and 2s ^S-2p 
^Pi/2 (1038 A) are predicted to differ by a factor of two 
(Aggarwal & Keenan 2004), then the line intensity ratio 
is: 

7(1032) ^ g,(1032)gfc(1038) 
7(1038) gfe(1032)g«(1038) 

where ga(k) is the probabihty for photon escape in 
the line-of-sight of the observer for line k and gb(k) is 
the angle-averaged escape probabihty. Both angle specific 
and angle -averaged escape probabilities (also known as 
escape factors) are required, with the former determining 



the transport of photons to the observer, which depends 
on the specific angle of observation. The angle-averaged 
escape factor controls the pumping of the upper state by 
photons from the rest of the plasma, which in turn deter- 
mines the number of excited states that can emit a photon, 
and this requires an average over all angles. It is the inclu- 
sion of both these factors that is needed to understand why 
the fine ratio can exceed the optically thin limit. 

At first sight one might think that if the upper state 
is pumped by the photon field, this involves the destruc- 
tion of the photon doing the pumping and would appear 
to suggest that there is no real gain in the population of 
the upper level. However, literature of the escape fac- 
tor method of accounting for photon absorption in a line 
or lines shows that this is not the case (see for example 
Mihalas 1970). Consequently, we maintain that our use 
of an escape factor approach coupled to a collisional ra- 
diative model takes into account all the relevant physics 
and is entirely realistic. This includes colUsional excita- 
tion (although de-excitation is much slower than spon- 
taneous radiative de-excitation as we are in collisional- 
radiative steady-state) and photo-excitation. Indeed, the 
latter is the process by which photons are absorbed by 
the plasma and that is accounted for by the escape factor 
method. 

For the case of a spherical plasma, ga(k) - gh(k) by 
symmetry, and hence 1(1032)/1(1038) = 2. However, for 
an infinite plane slab the factor becomes: 

7(1032) _ 
7(1038) ~ 

J^^d - exp[-^e^])dx f Ml - exp[-'fe^]}dndx 

f(l - exp[-fe-'])dx f //(I - exp[-^e^])dndx 

where tq is the optical depth of the 1038 A transition 
for a pathlength (Z) which is the perpendicular thickness 
of the plasma slab, and /i' = cos 6 where is the angle be- 
tween the perpendicular to the slab and the Une-of-sight to 
the observer (see Fig. 3). This expression is evaluated by 
numerical integration, and our results are shown in Fig. 
4 and discussed in Section 5. However, in Fig. 4 we plot 
the intensity ratio as a function of column density nj (i.e. 
product of electron density and pathlength) rather than tq 
as the former is a commonly-employed quantity and is 
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Fig. 1. Light curves for e Eri over the entire time period of the FUSE observations. The left-hand panel show the O vi 
light curve (1032 + 1038 A lines) and the right-hand one an expanded view of the first 25 ksec, which has the greatest 
level of stellar activity. Labels for phases 00 to 29, for intervals that had non-zero exposure of good signal-to-noise, are 
given and correspond to the Specnum names in the last column of Table 2; i.e. phase 02 did not have sufficient signal 
and is not in Table 2. The bin size is 200 sec for both Ught curves. Spectral bins are indicated by the dashed vertical 
lines. 
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Fig. 2. Light curves in the Ovi lines (1032 + 1038 A) for 11 Peg. The left-hand panel shows the entire time period of 
the FUSE observations, while the right shows an expanded portion of the individual observations, plus vertical dashed 
Unes for the exposure windows from IDF_CUT. Labels for phases pOO to p99 are given and correspond to the Specnum 
names in Table 2. The bin size is 200 sec for both light curves. 
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Fig. 3. Schematic diagram of the infinite plane slab ge- 
ometry discussed in Section 3. The line-of-sight to the 
observer is at an angle 6 to the perpendicular to the slab, 
which is of thickness /. 
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Fig. 6. Sample FUSE 1032/1038 line ratios in O vi plotted 
as a function of orbital phase. Spectra are in time order as 
compared to the first part of the eEri light curve and the 
total II Peg light curves. 



relatively straightforward to at least approximately esti- 
mate (see Section 4.4). 

Although our results formally only hold for an infi- 
nite slab, in reality they will be applicable to any case 
where the lateral dimension of the emitting plasma is 
much greater than the thickness /. Our theoretical models 
should therefore be a good approximation for the stellar 
transition regions considered here, given the small values 
of / found for these (see Sections 1 and 4.4), especially 
compared to the radius of a star. 

As the purpose of the current work is to illustrate the 
principle of this type of analysis we do not calculate the 
line intensity ratio for other geometries. However, we note 
that this will be the subject of a later paper. 

4. Results 

We have searched for line intensity enhancements in sev- 
eral active late-type stars observed with FUSE, employ- 
ing the the Ovi 1032/1038 line ratio as a diagnostic. 
Measured Ovi line fluxes and 1032/1038 intensity ra- 
tios, with their associated Icr errors, are listed in Table 
2. Below we discuss the measurements for each star and 
their variation as a function of orbital phase. 



4.1. eEri 

FUSE spectra for eEri extracted from the initial part 

of the observations showed 1032/1038 intensity ratios 
greater than the optically thin value (2.0) at statistically 
significant levels. For example, we find a measured line 
ratio of 2.58 ± 0.20 for Specnum 08 (see Table 2), which 



is larger than the theoretical value by 2.9cr, indicating a 
99.6% confidence level in our result being a true enhance- 
ment. We computed the orbital phase of each spectrum us- 
ing the ephemeris of eEri and rotational period of 11.35 
days from Croll et al. (2006), and the FUSE observations 
covered orbital phases from 0.90-0.93. The enhancement 
was observed for phases 0.903 and 0.93, and marginally 
for phase 0.909. However, spectra from intervals at later 
times (after 25 ksec) revealed no 1032/1038 ratios greater 
than 2.0. We show sample Ovi spectra for eEri in Fig. 
5, while the 1032/1038 ratios are plotted as a function of 
orbital phase in Fig. 6. 

4.2. II Peg 

II Peg showed several spectra with measured Ovi 
1032/1038 ratios greater than 2.0, but only a few are sta- 
tistically significant. For example, Specnum 04 has a mea- 
sured ratio of 2.30 ± 0.11 (see Table 2), larger than the 
theoretical value by 2.7cr, indicating a 99.3% confidence 
level in our result being a true enhancement. The Ovi 
light curve is nearly constant over the 300 ksec time in- 
terval of the observations (except for a small brightening 
in the middle set of data; 135-155 ksec), and the count 
rate in the FUSE band indicates the star was fairly ac- 
tive during this period. We computed the orbital phase of 
each spectrum using the ephemeris and rotational period 
of 6.72 days from Berdyugina et al. (1998). Our FUSE 
observations covered about 57% of the orbital period, in- 
cluding phases 0.44-0.50, 0.68-0.70 and 0.88-0.90. The 
1032/1038 intensity ratio was enhanced for phases be- 
tween 0.44-0.50, increasing from 2.14 to 2.30, and had 
a marginal enhancement at phase 0.91. We show sample 
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spectra for the O vi lines in Fig. 7, and the 1032/1038 ra- 
tios are plotted as a function of orbital phase in Fig. 6. 

4.3. Prox Cen 

Prox Cen (GJ 55 IC) was observed by FUSE in 2003 with 
the LWRS aperture. Christian et al. (2004) found several 
1032/1038 line ratios in the quiescent spectra which in- 
dicated the presence of significant opacity (i.e. values of 
< 2.0), and used these to derive pathlengths for the TR 
of up to ~ 13 km. However, these authors also measured 
1032/1038 ratios in two flare spectra which were >2.0, 
even allowing for observational uncertainities. These re- 
sults are hence also included in Table 2. 



4.4. Densities and patiniengtins 

To determine the column density of the Ovi region in 
our stellar sources, required for comparison with the- 
ory (see Section 3), we need estimates of both the den- 
sity and pathlength. Unfortunately, there are few elec- 
tron density {tie) diagnostics in the FUSE spectral re- 
gions for ions formed at similar temperatures to Ovi, 
which has a temperature of maximum fractional abun- 
dance in ionisation equilibrium of T^ - 10^'^ K (Bryans 
et al. 2009). The only possibility is Ne vi (formed at T^ = 



10^^ K; Bryans et al.), which has several density sensitive 
emission lines in the wavelength range -997-1010 A. 
However, the features are weak and noisy in the FUSE 
datasets for eEri and II Peg, and are not detected for Prox 
Cen. As a consequence, the Nevi lines provide at best 
only a constraint on the density range. For example, us- 
ing the latest version (V6.0) of the chianti database (Dere 
et al. 1997, 2009), the measured Ne vi 997.03/1005.7 ratio 
for II Peg of 1 . 1 +0.6 indicates an electron density in the 
range ~ 10"-10'^^ cm"^, while for eEri the experimental 
997.03/1005.7 ratio of 0.44+0.29 implies < lO'^cm^l 
These are consistent with the earlier results of Byrne et al. 
(1987) and Jordan et al. (2001), who derived - 10"^ 
and 10"^ cm"^ for II Peg and eEri, respectively, from a 
number of diagnostics. For Prox Cen, Giidel et al. (2002) 
found a typical electron density of ~ 10" cm"^ from O vii 
lines. Given these results, we adopt = 10" cm"^ for 
all 3 stars, which should be reliable to within an order of 
magnitude. 

Christian et al. (2006) have summarised estimates of 
the pathlengths (Z) for the O vi emitting regions in several 
late-type active stars, including Prox Cen. They find val- 
ues ranging from around I - 10-10^ km, and we therefore 
adopt a typical pathlength of I - 30 km, which should be 
accurate to an order of magnitude, similar to the electron 
density. Combining these, our expected column density 
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Phase Time O vi flux Ratio Comment 

MJD T„p (1032 A) (1038 A) 

(day) (sec) 10'" erg cm'^ s~' 10'" erg cm'^ s~' 1032/1038 Specnum 

eEri 

52993.4890 51816 0.160±0.003 0.080±0.003 1.93±0.08 Quiet 



0.903 


52992.7305 


867 


0.91±0.03 


0.54±0.03 


1.67+0.10 


00 


0.903 


52992.7383 


866 


1.44+0.03 


0.64+0.04 


2.26+0.16 


01 


0.909 


52992.8086 


659 


1.06±0.04 


0.41+0.03 


2.58+0.20 


08 


0.921 


52992.9375 


418 


1.02±0.04 


0.56+0.04 


1.81+0.15 


21 


0.922 


52992.9492 


801 


0.99±0.04 


0.60+0.04 


1.67+0.13 


22 


0.926 


52993.0001 


431 


0.67±0.03 


0.25+0.05 


2.72+0.54 


27, low S, 


0.927 


52993.0117 


866 


0.72±0.03 


0.37+0.03 


1.95+0.20 


28 


U.ylo 


jZyyi.y)lyj 


o2o 


1 IM\ 1 A AO 


A /I O 1 A AO 


z.z/±U.lj 


OA 

ly 


II Peg 
















51871.670 


15739 


3.03±0.04 


1.49+0.04 


2.03+0.07 


Total 


0.443 


51871.681 


1886 


2.92±0.08 


1.36+0.08 


2.14+0.14 


00 


0.454 


51871.750 


1919 


2.70±0.07 


1.22+0.07 


2.21+0.13 


02 


0.464 


51871.820 


1532 


2.68±0.05 


1.16+0.05 


2.30+0.11 


04 


0.680 


51873.272 


1074 


3.54±0.13 


1.68+0.18 


2.11+0.24 


49 


0.682 


51873.283 


786 


3.40±0.11 


1.42+0.12 


2.40+0.21 


50 


0.690 


51873.339 


687 


3.19±0.11 


1.54+0.12 


2.07+0.17 


51 


0.692 


51873.350 


1221 


3.48±0.11 


1.82+0.11 


1.91+0.13 


52 


0.701 


51873.416 


1496 


3.15±0.08 


1.50+0.09 


2.10+0.14 


54 


0.897 


51874.735 


1585 


2.83±0.08 


1.41+0.08 


2.00+0.12 


95 


0.908 


51874.804 


1805 


2.70±0.07 


1.28+0.08 


2.11+0.14 


97 


0.918 


51874.874 


1380 


2.57±0.09 


1.48+0.10 


1.74+0.14 


99 


Prox Cen 














Quiet 


52734.4748 


39840 


1.21±0.03 


0.66+0.03 


1.85+0.11 


Christian et al. 


Flare 1 


52734.4748 


600 


5.50±0.15 


2.44+0.15 


2.25+0.15 


Christian et al. 


Flare 2 


52735.1594 


1200 


0.97+0.05 


0.41+0.06 


2.40+0.35 


Christian et al. 



(2004) 
(2004) 
(2004) 



rigl (in units of cm ^) for the O vi plasmas in our stellar 
sample should be around log n^/ = 17.5. 



5. Discussion and future worl( 

First, we note that the simple detection of line intensity 
enhancements in several of our observations rules out a 
spherical geometry for the emitting plasmas, as theory 
predicts that a spherical region will show neither line in- 
tensity enhancement nor reduction, irrespective of the op- 
tical depth of the plasma (see Kerr et al. 2005 and Section 
3). This result assumes that the observer is not close to 
the plasma surface, which is of course a valid assumption 
for distant astronomical sources. Kerr et al. point out that 
their predictions for a spherical plasma may explain the 
few detections of opacity effects in stellar coronal emis- 
sion Unes, such as those of Fe xvii (Matranga et al. 2005), 
O vm and Ne x (Testa et al. 2007). One might expect high 
temperature coronal emission to be distributed (approx- 
imately) uniformly over the stellar surface, particularly 
under quiescent conditions, hence leading to a spherical 
geometry. 



To further illustrate how the detection of intensity en- 
hancement can constrain geometry we consider the the- 
oretical case outhned in Section 3, i.e. that of an infi- 
nite plane slab with column density thickness n^/, ob- 
served at an angle to the perpendicular. Fig. 4 shows 
the 1032/1038 intensity ratio of Ovi plotted as a func- 
tion of logarithmic column density for different values of 
6. From the figure, we see that significant intensity en- 
hancement is only predicted for a relatively narrow range 
of column density, log nj ^ 16.5-18.0, compatible with 
the approximate value derived for the stellar sources in the 
present paper, log n^/ ^ 17.5 (Section 4.4). Furthermore, 
the enhancement only occurs for small values of 9, indi- 
cating that the O vi emission regions are being observed 
close to face-on. In this regard, it is interesting to spec- 
ulate that the increase in the 1032/1038 ratio with time 
measured during the 0.4-0.5 orbital period of 11 Peg (Fig. 
6) may be due to the emitting plasma changing its ori- 
entation with respect to the line-of-sight such that 6 de- 
creases, resulting in an increase in the line intensity en- 
hancement. However, alternatively the plasma parameters 
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Fig. 5. Sample FUSE spectra of e Eri containing the 1032 and 1038 A transitions of O vi. The left-hand panel shows 
the quiet spectrum (1032/1038 intensity ratio = 1.93 + 0.08), while the flare spectra at several orbital phases are given 
in the right-hand one. Phases are ofi'set by 8x10"'^ erg cm"^ s"' for ease of presentation. 



may be changing, hence leading to an increase or decrease 
in column density. 

It can therefore be seen that the observation of line in- 
tensity enhancement in astronomical spectra due to opac- 
ity does allow some constraints to be placed on both the 
geometry of the emitting plasma and orientation to the ob- 
server. However, further work is required. For example, it 
would be interesting to properly observe a stellar source at 
high spectral and temporal resolution over a full rotation 
period, to investigate if the line intensity enhancement can 
be coiTelated with rotation. This would provide evidence 
that the emitting region is changing its orientation with 
respect to the sightline. It would also be useful to extend 
the work to astronomical sources other than cool stars. 

Finally, although our research provides evidence for 
line intensity enhancements due to opacity in astrophys- 
ical sources, it would be highly desirable to confirm the 
theory of Kerr et al. (2004, 2005) in a well-diagnosed lab- 
oratory plasma. In particular, using a laboratory plasma 
of known geometry will allow a strigent test of the theory 
in terms of plasma geometry and orientation of the ob- 
server, as different results are expected depending on e.g. 
whether the plasma is slab, cylindrical or spherical (Kerr 
et al. 2004, 2005). 

Here we provide an outline design of an experiment 
using a high-power laser to check our theoretical mod- 
elling, through the observation of line ratios from a lab- 
oratory plasma that has been independently diagnosed. 
The experiment proposed below builds on recent work 



conducted using the NIKE laser (Back et al. 2006). In 
these studies. Back et al. demonstrated that it was possible 
to heat a low density (few mg cm"-') aerogel foam target 
to electron temperatures of over 5x10^ K with good uni- 
formity. This was accomplished by driving a supersonic 
heating wave through the foam, producing a plasma of 
known density that disassembles through the progress of 
a rarefaction wave starting at the outside of the plasma. 
By doping the foam with an element in a central region 
of chosen geometry, a plasma of known, uniform tem- 
perature (diagnosed using Thomson scattering), uniform 
initial density, known composition and selected geometry 
can potentially be produced in the laboratory. The obser- 
vation of a line intensity ratio (optically thick/optically 
thin) as a function of increasing dopant concentration and 
angle will allow us to produce in the laboratory the ana- 
logue of the astrophysical situation described in KeiT et 
al. (2004, 2005), and thereby check our modelling. 

The schematic diagram of the experiment is shown in 
Fig. 8. A laser heats a block of aerogel foam at sub-critical 
density into which is doped the emitting material (for this 
excercise we choose the doped material to be sodium) 
in a region of a chosen (in this case planar) geometry. 
Thomson scattering determines the electron and ion tem- 
perature (both are needed: the first determines the kinetics 
and the second determines the line width and thereby the 
optical depth). The doped region remains unaffected until 
it is disrupted by the rarefaction wave that moves inwards 
from the outside of the foam. Calculations with the NIMP 
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Fig. 7. Sample FUSE spectra of II Peg containing the 1032 and 1038 A transitions of O vi. The left-hand panel shows 
the total II Peg FUSE spectrum (1032/1038 intensity ratio = 2.03 + 0.07), while selected orbital phases are given on 
the right. Phases are offset by 6x10"'-' erg cm"^ s"' for ease of presentation. 
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Fig. 8. Schematic diagram of a possible experiment to 
measure line intensity enhancements arising from opacity 
effects in a laboratory plasma. See Section 5 for further 
details. 



code of Rose (1997) indicate that for sodium doped into 
the central region, the kinetics reach a steady state in a 
few hundreds of picoseconds, which shows that the doped 
plasma has sufficient time to reach a steady-state before it 
is disrupted. The spectroscopic observation of the line in- 
tensity ratio is conducted at a variety of angles (0) and for 
a number of different dopant concentrations. 



Acknowledgements. FPK is grateful to AWE Aldermaston for the award 
of a William Penney Fellowship. DC thanks the Royal Society for an 
International Travel Grant in support of this research, chianti is a collab- 
orative project involving researchers at the Naval Research Laboratory 
(USA), Rutherford Appleton Laboratory (UK), and the Universities of 
Cambridge (UK), George Mason (USA) and Florence (Italy). 



References 

Aggarwal, K. M., & Keenan, F P. 2004, Phys. Scripta, 70, 222 

Back, C. A., Feldman, U., Weaver, J. L., et al. 2006, 

J. Quant. Spec. Radiat. Transf , 99, 21 
Berdyugina, S. V., Jankov, S., Ilyin, I., Tuominen, I., & Fekel, F. C. 

1998, A&A, 334, 863 
Bloomfield, D. S., Mathioudakis, M., Christian, D. J., Keenan, F P, & 

Linsky, J. L. 2002, A&A, 390, 219 
Bryans, R, Landi E., & Savin, D. W. 2009, ApJ, 691, 1540 
Byrne, R B., Doyle, J. G., Brown, A., Linsky, J. L., & Rodono, M. 1987, 

A&A, 180, 172 

Christian, D. J., Mathioudakis, M., Bloomfield, D. S., Dupuis, J., & 

Keenan, F P 2004, ApJ, 612, 1140 
Christian D. J., Mathioudakis, M., Bloomfield, D. S., Dupuis, J., 

Keenan, F P, Pollacco, D. L., & Malina, R. F 2006, A&A, 454, 

889 

CroU, B., Walker, G. A. H., Kuschnig, R., et al. 2006, ApJ, 648, 607 
Dere, K. P., Landi, E., Mason, H. E., Monsignori-Fossi, B. C, & Young, 

R R. 1997, A&AS, 125, 149 
Dere, K. R, Landi, E., Young, P R., Del Zanna, G., Landini, M., & 

Mason, H. E. 2009, A&A, 498, 915 
Dixon, W. v., Sahnow, D. J., Ban-ett, P E., et al. 2007, PASR 1 19, 527 
Gillaspy, J. D., Lin, T., Tedesco, L., Tan, J. N. et al. 201 1, ApJ, 728, 132 
Gudel, M., Audard, M., Skinner, S. L., & Horvath, M. I. 2002, ApJ, 580, 

L73 

Jordan, C, Sim, S. A., McMurry, A. D., & Aruvel, M. 2001, MNRAS, 
326, 303 



10 



F. P. Keenan et al.: Intensity enhancement due to opacity 



Kerr, F. M., Rose, S. J., Wark, J. S., & Keenan, F. P. 2004, ApJ, 613, 
L181 

Kerr, F M., Rose, S. J., Wark, J. S., & Keenan, F P. 2005, ApJ, 629, 
1091 

Ko, Y.-K., Raymond, J. C, Zurbuchen, T. H., Riley, P., Raines, J. M., & 

Strachan, L. 2006, ApJ, 646, 1275 
Kohl, J. L., & Withbroe, G. L. 1982, ApJ, 256, 263 
Mathioudakis, M., McKenny, J., Keenan, F. P., WilUams, D. R., & 

PhilUps, K. J. H. 1999, A&A, 351, L23 
Matranga, M., Mathioudakis, M., Kay, H. R. M., & Keenan, F. P. 2005, 

ApJ, 621, L125 

MihaUs, D. 1970, Stellar Atmospheres, W. H. Freeman & Co. Ltd. 
Moos, H. W., Cash, W. C, Cowie, L. L., et al. 2000, ApJ, 538, LI 
Nakagawa, A. 2008, ApJ, 674, 1167 

Rose, S. J., 1997, Rutherford Appleton Laboratory Technical Report 
RAL-TR-97-020 

Rose, S. J., Matranga, M., Mathioudakis, M., Keenan, F. P., & Wark, J. 

S. 2008, A&A, 483, 887 
Sahnow, D. J., Moos, H. W., Ake, T. B., et al. 2000a, ApJ, 538, L7 
Sahnow, D. J., Moos, H. W, Ake, T. B., et al. 2000b, Proc. SPIE, 4013, 

334 

Testa, P, Drake, J. J., Peres, G. , & Huenemoerder, D. P 2007, ApJ, 665, 
1349 



